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Abstract. We formulate an extended linear a model of a quarkonia nonet and a tetraquark nonet as 
well as a complex iso-singlet (glueball) field to study the low-lying scalar meson. Chiral symmetry 
and Ua(1) symmetry and their breaking play important role to shape the scalar meson spectrum in 
our work. Based on our study we will comment on what may be the mass of the lowest possible 
scalar and pseudoscalar glueball states. We will also discuss on what may be the nature of the sigma 
or /o(600) meson. 
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INTRODUCTION 



Theoretical understanding of the so-called Higgs Boson of QCD namely the a or 
/o(600) meson is important as they play an important role in chiral symmetry break- 
ing and can be used as a probe for QCD vacuum. While the confirmation of its existence 
from %% scattering process puts an end to the decades-long controversy [1,2], the agree- 
ment on the nature of a//o(600) however has not yet been achieved. The tetraquark 
model [3] supports cr//o(600) to be a tetraquark state. Whereas, recent data from %n 
and 77 scattering [4, 5] suggest it has a sizable fraction of glueball and another study 
from K-matrix analysis [6] suggested that /o(600)/<7 should be a glueball dominant 
state. 

On the other hand the agreement on the lowest possible iso-scalar and pseudoscalar 
glueball states has neither been acheived. The result from Lattice simulations suggests 
/o(1500) or /o( 1700) could be a glueball rich iso-scalar, while T](1489) can be a glueball 
rich pseudoscalar [7]. 

The great success of chiral symmetry in understanding the nature of lightest pseudo- 
scalars motivates us to use an extended version of linear o models [8, 9] to study the 
nature of these low-lying isoscalars and pseudoscalars. The motivation for constructing 
an extended linear sigma model consisting of effective quarkonia, tetraquark and glue- 
ball fields comes from the physical considerations that scalar condensates are allowed 
by the QCD vacuum. So in principle apart from quarkonia and tetraquark condensates, 
scalar glueball condensate should also be present in the QCD vacuum and need to be 
considered while studying the vacuum excitation. The mixing pattern followed in our 
framework comes from the consideration that mesons having identical external quan- 



turn numbers can mix even if they have different internal flavour structures. Based on 
these premises we attempt to study the nature of cr// (600) meson as well as lowest iso- 
scalar and pseudoscalar glueball candidate. In the following section we briefly review 
the model Lagrangian and the symmetry breaking pattern. At last we will present our 
results and conclude. 

MODEL 

The extended linear o model can be systematically formulated under the symmetry 
group SUr(3) x SUl(3) x Ua(1). Three types of chiral fields are included: a 3 x 3 matrix 
field $ which denotes the quarkonia states, a 3 x 3 matrix field <J»' which denotes the 
tetraquark state, and a complex chiral singlet field Y which denotes the pure glueball 
states. 

Up to the mass dimension 0(p A ) (we assume that they are the most important op- 
erators to determine the nature of light scalars of ground states), the Lagrangian of 
our model can include two parts: the symmetry invariant part and the symmetry 
breaking one J^sb- The symmetry invariant part includes those terms which respect 
SU(3)l x SU(3)r symmetry as well as Ua(1) symmetry: 

£> s = Tr(d^d^) + Tr(d^'d^') + d^Yd^Y* - m<i , 2 Tr{&<$>) - m^ 2 Tr{^'^') 
-m Y 2 YY*- AiTr(4) t 4>4> t 4>) - XiTr(&'<l>'&'<l>') - X 2 Tr - X Y (YY*) 2 

-[X 3 £ abc £ de f<S> d a <S> e b <S>/ c + h.c.} + [kYDet(<$>) +h.c], (1) 

While the symmetry breaking part includes the following terms 

J^ss = [7>(fl.4>) +h.c] + [Tr(B'.<S>') +h.c] + (D.Y + h.c.) - [X m Tr(<i><t> p ) +h.c.]. (2) 

To construct the symmetric terms in the Lagrangian we closely follow [10] where the 
choice of terms are limited by the number of internal quark plus antiquark lines at the 
effective vertex, which is set to 8. This condition is relaxed for two terms with coupling 
constants X\ and X2 and the reason to include these two terms stems from the practical 
consideration of making sure our potential is bounded. For detailed discussion on the 
terms appearing in the Lagrangian please refer to [11]. As evident from the Lagrangian 
both chiral symmetry SU (3)l x SU (3)r and U(\)a symmetry are explicitly broken by 
the terms in S£sb- The 3x3 matrices B and B' responsible for the breaking of the 
symmetry can be parametrized as: B(B') = T a b a (T a b a ') (T a are the generators of U(3) 

with T = 

Since the vacuum expectation values of the quarkonia and tetraquark fields can carry 
those quantum numbers which are allowed by the QCD vacuum, only a = (0, 3, 8) fields 
are allowed. The choice of this external fields b a (b a ') control the nature and extent of the 
symmetry breaking. Out of various possible symmetry breaking scenarios, we consider 
the following case in our present study where: 

. b Q (b ') ^ 0, b 3 (b 3 ') = and b$(b$') ^ 0. In this case, U(3) A and SU(3) V both 
are explicitly broken and SU(3)v is broken to SU(2)y. As a result m u = m c i^m s , 



TABLE 1. The categorization of scalar and pseudo-scalar states in term of 
isospin quantum number are demonstrated. States in the same category can mix 
with each other. 



Isospin 


7=1 




7 = 


PseudoScalars(P=- 1 ) 




{K,K'}, {K*,K*'} 


{T?1,T?2,T?3,7?4,TJ 5 } 


Scalars(P=l) 


{a, a 1 }, 


{k, k'}, {k*, k*'}, 


{/l!/2,/3)/4,/5} 



where m, is the quark mass of the i th flavour. 

This is reasonable considering the up and down quark masses are nearly equal to each 
other and thereby indicating SU (2)y is a good (approximate) symmetry. The remnant 
SU(2)y isospin symmetry allows us to represent two condensates for quarkonia and 
teraquark fields each as: vo, vg and v'o, v'% respectively. While the gluonic condensate in 
our theory is labeled as v y . 

RESULTS AND CONCLUSION 

Due to the unbroken SUy(2) isospin symmetry, physical scalar and pseudo-scalar states 
can be categorized into three groups with isospin quantum numbers as 7" = 1 (triplet), \ 
(doublet) and 0, respectively. Only bare quarkonia, tetraquark and glueball fields with 
the same isospin quantum number can mix with each other to form physical states. 
Moreover, there is no mixing between scalar and pseudoscalar fields. Thus the chiral 
singlet glueball field can only mix with the isospin singlets of quarkonia and tetraquark 
fields. Using these facts, the physical states below 2 GeV can be tabulated as given in 
Table 1, where the isodoublet {K,K'} is connected with the isodoublet {K*,K*'} by 
charge conjugation. And a similar relation holds for {k, k'} and {k*, k*'}. 

There are 15 parameters in our model. To solve these parameters we treat the 
tetraquark vacuum condensates as well as the mixing angles for isotriplet (9 K ) and 
isodoublet (6k) fields as input parameters. Then from the mass matrices of a, a', K and 
k', parameters related to isotriplet and isodublet sectors are solved. The symmetry pa- 
rameters {bo, bg, b' Q , b'g,D} are solved from the vacuum stability conditions. Whereas, to 
solve for the parameters related to glueball sector, we use following two conditions: 

Model — 

Tr[M n 2 ] Exp , Det[M 2 \ 

Model 

= Det[M n 2 ] Exp . (3) 

At the end we are left with one free parameter which is bare glueball mass my and is 
used in our study as a scanning parameter. 

For space constriant we refer to [11] for details on parameter fixing and method 
to choose the best fit solution. In choosing the best fit solution we vary the n' mass 
from 1.2—1.4 GeV and choose those solution which give the tree level decay width for 
/o(600) ->• n% between 0.35 - 0.9 GeV. 

Our best fit parameter set is presented in Table (2). We would like to highlight a 
few features out from it. 1) In absence of the explicit symmetry breaking terms, it is 
the negative mass parameter m$ 2 that would trigger the spontaneous chiral symmetry 
breaking. 2) The sign of vy is correlated with the sign of k, and the sign of k is determined 
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FIGURE 1. (Left figure) The dependence of Xy upon m Y is demonstrated. A solid circle marker shows 
the point Xy = 0, which corresponds to m Y — 3.452 (my — 1 .858) . (Right figure) The dependence of mass 
of fo upon my is demonstrated. A vertical line with m Y is drawn to read out the lowest mass m f s = 1 .86. 



TABLE 2. The values of parameters in our fit are shown where the best 
value of m K i is found to be m K i = 1 .2 GeV. 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 


9 n (radian) 


-0.604 


V 


8.248 


m® 2 {GeV 2 ) 


-0.025 


9k (radian) 


-0.714 




76.428 


b ' {GeV 3 ) 


0.166 


vo (GeV) 


0.074 


X 3 (GeV) 


-0.738 


m<j>/ 2 {GeV 2 ) 


0.744 


v 8 (GeV) 


-0.115 


Ay 


38.327 


h' {GeV 3 ) 


0.18 


v ' (GeV) 


0.203 


k 


-78.15 


h 


35.465 


v 8 ' (GeV) 


0.126 


Kn {GeV 2 ) 


-1.044 


D {GeV 3 ) 


-0.265 


v y (GeV) 


-0.109 


bo {GeV 3 ) 


-0.085 






my 2 {GeV 2 ) 


3.0 


b s {GeV 3 ) 


-0.161 







from the mass spectra of pseudoscalar sector. 3) The couplings X\, X[, Xy are positive 
which guarantee the potential is bounded from below. 4) The values of X\, Xj, and Xy as 
well as k are large, which demonstrate the non-perturbative nature of the model. 

It is found that the condition Xy > can predict the lightest glueball scalar should 
be around 2.0 GeV or so, as can be read off from Fig. (lb), while the lightest glueball 
pseudo scalar should be 775. The mass splitting between these two glueball states is 
controlled by parameters vy and Xy and is found to be around 0. 15 GeV. When compared 
with the Lattice QCD prediction for the glueball bare mass reported in [12] where the 
mass is 1.611 GeV, our result my = 1.73 GeV is slightly heavier than this prediction. 
When my = 1.611 GeV is taken, then the predicted mass of the lightest glueball is 
mjo = 2.29 GeV. 3) The lightest scalar / 1 °(600) is found to be 0.27 GeV or so and is a 
quarkonia dominant state. 

In Figure 1, we demonstrate the dependence of Xy and f° masses upon the free 
parameter m y with the rest of parameters are given in Table (2). As shown in Fig. 
(la), when my is larger than 3.4 GeV 2 , the Xy becomes negative. Then the potential 
of our model has to confront with the problem of unbounded vacuum from below. In the 
allowed values of m\, the masses of jf , i = 1,2,3,4 are almost independent of its value, 



TABLE 3. Mass spectra and components for the pseudo-scalar amd scalar mesons based on our 
fit are shown where the best value of m n i is found to be m n i = 1.2 GeV. 



Meson 


Our Value 

(GeV) 


quarkonia (%) 


tetraquark (%) 


glueball (%) 


Experimental 
Value (GeV) 


'15 


1 858 


037 


001 


99 96? 


1 756 + 009 


14 


1.380 


75.803 


24.167 


0.03 


1.476 ±0.004 




1.291 


26.700 


73.294 


0.006 


1.294 ±0.004 


m 


0.907 


15.852 


84.145 


0.003 


0.958 ± 24 x 10~ 5 


m 


0.595 


81.607 


18.393 


0.0 


0.548 ± 24 x 10~ 6 


fs° 


2.09 


0.01 


0.0 


99.99 




/4° 


1.487 


77.469 


22.53 


0.001 


1.505 ±0.006 


/3° 


1.347 


22.177 


77.82 


0.003 


1.2-1.5 


/2° 


1.124 


21.561 


78.439 


0.0 


0.980 ± 0.010 


/l° 


0.274 


78.784 


21.211 


0.005 


0.4-1.2 



as demonstrated in Fig. (lb). The upper bound of niy is determined from the condition 
r /0 > 0.35 GeV. 

To conclude, our model predicts that the isoscalar glueball should be heavier than 
2.0 GeV when the pseudoscalar 77(1726) is the best glueball candidates. The lowest 
isoscalar /o(600) is found to be quarkonia dominant state with a considerable tetraquark 
component. 
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